arXiv: 1504.02108v2 [cond-mat.str-el] 11 Sep 2015 


Scattering of high-energy magnons off a magnetic skyrmion 

Sarah Schroeter and Markus Garst 

Institut fur Theoretische Physik, Universitat zu Koln, Zulpicher Str. 11a, 50931 Koln, Germany 

(Dated: September 14, 2015) 

We discuss the scattering of high-energy magnons off a single magnetic skyrmion within the field- 
polarized ground state of a two-dimensional chiral magnet. For wavevectors larger than the inverse 
skyrmion radius, kr s 1 , the magnon scattering is dominated by an emerging magnetic field whose 
flux density is essentially determined by the topological charge density of the skyrmion texture. This 
leads to skew and rainbow scattering characterized by an asymmetric and oscillating differential cross 
section. We demonstrate that the transversal momentum transfer to the skyrmion is universal due to 
the quantization of the total emerging flux while the longitudinal momentum transfer is negligible in 
the high-energy limit. This results in a magnon-driven skyrmion motion approximately antiparallel 
to the incoming magnon current and a universal relation between current and skyrmion-velocity. 

PACS numbers: 


I. INTRODUCTION 

The experimental discovery of skyrmions in chiral 
magnet 7 and in magnetic monolayer^ 10 has trig¬ 
gered an increasing interest in the interaction of spin 
currents with topological magnetic texturesJUMU p. p as 
been demonstrated 13 16 that skyrmions can be manip¬ 
ulated by ultralow electronic current densities of 10 6 
A/m 2 , which is five orders of magnitudes smaller than in 
conventional spintronic applications using domain walls. 
The adiabatic spin-alignment of electrons moving across 
a skyrmion texture results in an e mergent electrodynam¬ 
ics implying a topologica l 11 * 12 * 3 ^ as well as a skyrmion- 
flow Hall effect!- 17 In insulators, the interplay of thermal 
magnon currents and skyrmions is marked by a topolog¬ 
ical magnon Hall effect and a magnon-driven skyrmion 
motion. 23 -^* The topological nature of the magnetic 
skyrmions is responsible for a peculiar dynamics 33 37 1 
that is also at the origin of these novel spintronic and 
caloritronic phenomena, which are at the focus of the 
fledgling field of skyrmionicsW^ 

In two spatial dimensions, skyrmions are identified by 



FIG. 1: (a) A chiral magnetic skyrmion texture of linear size 
r s . (b) Illustration of a classical magnon trajectory within 
the x-y plane scattering off a skyrmion positioned at R with 
impact parameter b and classical deflection angle 0. 


the topological charge density 

Ptop = 2-n(a x n x dyh), (1) 

where n is the orientation of the magnetization vec¬ 
tor. For a magnetization homogeneously polarized at the 
boundary, the spatial integral f d 2 rp top = W is quan¬ 
tized, W G Z, and thus allows to count skyrmions within 
the sample. In turn, a finite winding number W trans¬ 
lates to a gyrocoupling vector G in the Thiele equation 
of motion of the skyrmionf 38 and the resulting gyrotropic 
spin-Magnus force governs its dynamics.^ As a conse¬ 
quence, in the presence of an applied electronic spin cur¬ 
rent, the skyrmions will acquire a velocit y 4 * 15 * 17 * that re¬ 
mains finite in the limit of adiabatic spin-transfer torques 
and small Gilbert damping a , giving rise to a universal 
current-velocity relation.^* 

In order to address the interaction of magnon cur¬ 
rents with magnetic textures, a corresponding adiabatic 
approximation has been recently invoked on the level 
of the Landau-Lifshitz-Gilbert equation by Kovalev and 
TserkovnyakPSl This approximation has been used in 
Refs. 1231241 to derive an effective Thiele equation of mo¬ 
tion for the skyrmion coordinate R in the presence of a 
magnon current density J, 

G x R = — G x + f3v e ff + ..., (2) 

with /3 = 0 in the adiabatic limit. The effective velocity 
u eff = gp B J/(hm 0 ) is related to the current density via 
the g-factor g , the Bohr magneton (jlb > 0 and the lo¬ 
cal magnetization m q. The gyrocoupling vector is given 
by G = —ATTzhmo/^ggB) with units of spin density cor¬ 
responding to a flux of —27 xh per area of a spin-^ in 
a two-dimensional system with the unit normal vector 
z. The dots in Eq. © represent further terms omitted 
for the purpose of the following discussion, that is, in 
particular, a damping force proportional to the Gilbert 
constant a. Neglecting these additional terms, Eq. © 
predicts for /3 = 0, similar to the skyrmion-driven mo¬ 
tion by electronic currents, a universal current-velocity 
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relation R = — v e ^ = —giiBJ/{hm o) with a skyrmion ve¬ 
locity that is antiparallel to J. Consequently, a magnon 
current generated by a thermal gradient will induce a 
skyrmion motion towards the hot region of th e sample, 
which was indeed observed numericallyf 23 *^- 2 -^ Mochizuki 
et alW ^ also used Eq. © with P a= 0 to account for the 
experimental observation of a thermally induced rotation 
of a skyrmion crystal. 

However, the question arises as to when the adiabatic 
limit of Eq. (2| is actually applicable and under what 
conditions. The validity regime of the adiabatic approx¬ 
imation for magnon-driven motion of magnetic textures 
has not been explicitly discussed in Ref. [40j In fact, in 
order to account quantitatively for their numerical exper¬ 
iment Lin et al 24 introduced the ft parameter in Eq. © 
on phenomenological grounds calling it a measure for 
non-adiabaticity. Subsequently, Kovalev 2 ® argued that 
a finite (3 parameter arises due to dissipative processes. 

In contrast, we have recently shown by considering the 
magnon-skyrmion scattering problem 29 that a monochro¬ 
matic magnon current with energy 5 will give rise to a 
reactive momentum-transfer force in the Thiele equation 
which reads in linear response 


G X R — /ccr_i_ (k) (/£ x J £ ^j T /C(T|| (k) t/ £ T ..., (3) 

where the magnon dispersion is 5 = £ ga ,p + (hk) 2 / (2M mag ) 
with the magnon gap £ gap and the magnon mass M mag . 
This force on the right-hand side of Eq. © is determined 
by the two-dimensional transport scattering cross sec¬ 
tions 



1 — cos x 
-siny 


da 

dx 


( 4 ) 


where is the energy-dependent differential scattering 
cross section of the skyrmion. In the limit of low energies 
kr s < 1, where r s is the skyrmion radius, s-wave scat¬ 
tering is found to dominate so that a±(e) —>> 0 and, as 
shown in Ref. (29j the force becomes longitudinal to J £ . 
This, in turn, implies a skyrmion motion approximately 

perpendicular to the magnon current, R —)> z x X, 

\G\ 

thus maximally violating the predictions of the adiabatic 
limit of Eq. ©• This implies that Eq. Q is not valid for 
low-energy magnons whose wavevector is comparable or 
smaller than the inverse size of the texture. 

It is one of the aims of this work to demonstrate ex¬ 
plicitly that in the high-energy limit, kr s 1, on the 
other hand, the momentum-transfer force of Eq. © due 
to a monochromatic magnon wave indeed reduces to the 
form of Eq. ©• The effective velocity in this case, how¬ 
ever, is to be identified with v e ^ = \A\ 2 hk/M mag where 
A is the amplitude of the incoming magnon wave. In 
the high-energy limit the magnon-skyrmion interaction 
is dominated by a scattering vector potential, i.e., an 
emerging orbital magnetic field whose flux is quantized 
and related to the skyrmion topology. As a result, the 


transversal momentum transfer assumes a universal value 
in the high-energy limit ka±(e) An as anticipated in 
Ref. 25, Moreover, the longitudinal momentum trans¬ 
fer yields a reactive contribution, /3 e , to the /3 parameter 
that, in this limit, is determined by the square of the clas¬ 
sical deflection function 0(5) integrated over the impact 
parameter 6, see Fig. EH 

Pe = ^ k j™Jb{Q(b)) 2 . (5) 

As the scattering is in forward direction at high energies, 
0(6) 1/fc, the parameter vanishes as [3 £ oc 1/k so that 

it is indeed small for large kr s 1. 

The outline of the paper is as follows. In section |Tl 
we shortly review the definition of the magnon-skyrmion 
scattering problem and some of the main results of 
Ref. ESi In section |III| we examine the scattering prop¬ 
erties of high-energy magnons including the skew and 
rainbow effects, the total and transport scattering cross 
sections, and the magnon pressure on the skyrmion lead¬ 
ing to Eq. pi). We finish with a short discussion in section 

m 


II. SKYRMIONIC SOLITON AND ITS 
SPIN-WAVE EXCITATIONS 


This section closely follows Ref. |29| and reviews 
the magnon-skyrmion scattering problem in a two- 
dimensional chiral magnet. We start with the stan¬ 
dard model for a cubic chiral magnet restricted to a 
two-dime nsiona l plane that is described by the energy 
functiona l 41 * 42 * 


c _ Ps 

^ - Yl 


(daUj ) 2 


~ 2ft 2 flB 


(6) 


with spatial index a G {1, 2} = {x, y} and i, j G {1, 2, 3}, 
Ci a j is the totally antisymmetric tensor with 6123 = 1, 
and p s is the stiffness. The two length scales are given 
by the wavevectors Q and n. The former determines 
the strength of the spin-orbit Dzyaloshinskii-Moriya in¬ 
teraction, that we chose to be positive, Q > 0. The 
latter, n > 0, measures the strength of the applied mag¬ 
netic field, that is applied perpendicular to the two- 
dimensional plane, B = z. We neglect cubic anisotropies, 
dipolar interactions as well as magnetic anisotropies for 
simplicity. The latter can be easily included resulting in 
an additional length scale. 


A. Skyrmionic saddle-point solution 

The theory (J6| possesses a topological soliton solu¬ 
tion, i.e., a skyrm ion, as first pointed out by Bogdanov 
and Hubert With the standard parametrization of 
the unit vector hj = (sin#cos(^, sin^sin^, cos#), the 
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skyrmion obeys 


0 = 6(p), v? = x+|, (7) 


where p and x are polar coordinates of the two- 
dimensional spatial vector r = p(cosy, siny). The polar 
angle 9 obeys the differential equation 


9 H-■— 


9' sin 9 cos 9 2 Q sin 2 9 


— ft 2 sin# = 0 , 


(8) 


with the boundary conditions 0(0) = 7r and 
linip^oo 9(p) = 0. At large distances pK 1, the po¬ 
lar angle obeys the asymptotics 9(p) ~ e~ Kp / ^fp, which 
identifies ft as the inverse skyrmion radius. The resulting 
skyrmion texture is illustrated in Fig. [lja). The associ¬ 
ated topological charge density 

Ptop = ^n s (d x n 3 x d y n s ) = — (9) 

integrates to f d 2 rp^ op = — 1 identifying the solution as 
a skyrmion. The skyrmion radius r s can be defined with 
the help of the area f d 2 r(l — h z )/2 = nr 2 , and it is found 
to approximately obey r s ~ 1/ft 2 . 

The skyrmion is a large-amplitude excitation of the 
fully polarized ground state as long as its energy is pos¬ 
itive, which is the case for ft > ft cr where ft 2 r ~ 0.8Q 2 , 
which is the regime we focus on. For smaller values of 
ft, skyrmions proliferate resulting in the formation of a 
skyrmion crystal ground state. 


phase factor —e~ %x in the second term. For the discussion 
of magnon scattering, it will be convenient to introduce 
a wavefunction ipi^ with respect to a frame that reduces 
to the laboratory frame at large distances, that is simply 
obtained by the gauge transformation 

VWM) = -e~ lx ip(r,t). ( 12 ) 


2. Magnon Hamiltonian 

In order to derive an effective Hamiltonian for ip, we 
consider the Landau-Lifshitz equation 

d t n = -yn x H eff , (13) 

with 7 = gpB/h, where the effective magnetic field 
^eff(r,£) = si^t) determined by the functional 
derivative of the integrated energy density E = f dtdr£. 
Expanding ( jl3| ) in lowest order in ip, one finds that 
the spinor = (ip,ip*) is governed by a bosonic 
Bogoliubov-deGennes (BdG) equation 


ihr z d t $ = W&, 


(14) 


with the Hamiltonian 


H = 


h 2 (—ilV — r^a ) 2 

2 M ma g 


+ IVo + T X V X , 


(15) 


B. Magnon-skyrmion scattering problem 

1 . Magnon wavefunction 

The magnons correspond to spin-wave excitations 
around the skyrmion solution n s that can be ana¬ 
lyzed in the spiri t of previous work by Ivanov and 
collaborators.^^ We introduce the local orthogonal 
frame e^ej = Sij with q x e 2 = e 3 , where 63 ( 1 *) = 
n s ( r) tracks the skyrmion profile. For the two or¬ 
thogonal vectors we use ej = (— sin<^, cos p, 0) and 
e 2 = (— cos 9 cos p, sin 9 sin p, cos 9) . The excitations are 
parametrized in the standard fashion 

n = e^\/l - 2|^| 2 + e + ip + e_ip*, (10) 

where ip is the magnon wavefunction and e± = ~p^(&i db 
ie 2 ). For large distances, p r s , this parametrization 
assumes the form 


where V T = (d x ,d y ), and r x and r z are Pauli matrices. 
The potentials are given by 


V„W = SgE(- 


sin 2 9 Q sin(2 9) 


2 p 2 


2 P 




91 A v 

- Q 2 sin 2 9 + k 2 cos 9 — Q9’ —— J 
sin 2 9 Qsin(26») , 


9 12 ' 
2 

_ ^ 

2 /’ 


(16) 


(17) 


The magnon energy gap is defined by 


5 gap — 


gpBpsK 


m 0 


H 2 k 2 

2M ma o 


(18) 


which also identifies the magnon mass M mag . The vector 
potential reads a = a x (p)x with y T = (— siny, cosy) 
and 


a /1 - 2 |- 0| 2 + (-j={% + iy)(~e lx ip)+c.c.). 


V V2 


( 11 ) 


It is important to note that the local frame corresponds 
to a rotating frame even at large distances reflected in the 


a x = ( ^--Qsm9. (19) 

P 

It obeys the Coulomb gauge Va = 0. The polar angle 
in all potentials is the soliton solution, 9 = 9(p), and 
depends on the distance p. 
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FIG. 2: Regular part of the effective magnetic flux density 
(24) for various values of ft 2 /Q 2 . For lower values of k 2 /Q 2 
the flux density close to the skyrmion center is suppressed 
and even becomes negative for k 2 /Q 2 < 1.3. As a result, 
the effective local Lorentz force evaluated along a classical 
magnon trajectory with b — 0 changes sign resulting in a 
suppression of the deflection angle. 


3. Effective magnetic flux 


Far away from the skyrmion the Hamiltonian simplifies 
EL —ELo for p —>> oo with 


EL 0 = 




r zl a x ) 2 


2 M, 


+ le, 


gap- 


( 20 ) 


mag 


The remaining vector potential is attributed to the choice 
of the rotating orthogonal frame in the definition of the 
magnon wavefunction, see Eq. (11). It can be easily elim¬ 


inated by the gauge transformation (12), 


T \F 


lab 


«fab = a 


e ~iT z (x+n) iff 

„ 1 cos 0 

p 


p 


Qsin#. 


( 21 ) 

( 22 ) 


With respect to this laboratory orthogonal frame, the 
vector scattering potential ai a b = CL\ ah X vanishes expo¬ 
nentially for large distances, p r s . 

The associated flux B = V x (fiai a b) = Bz will play 
an important role in the following discussion, where 
B(r) = ^d p (paf ah (p)). According to Stokes’ theorem the 
total flux f d 2 rB(r) = 0 vanishes as ai a b is exponentially 
confined to the skyrmion radius. However, there is an 
interesting spatial flux distribution, 


B( r) = —47rM(r) + B reg (|r|), (23) 

Breg(p) = ^-Ptop - sin6 o) . (24) 

Since for small distances ^ ah (p) —> —2/p, there is a sin¬ 
gular flux contribution at the skyrmion origin with quan¬ 
tized strength —Anh. As it is quantized, this singular 
flux will not contribute to the magnon scattering. The 
regular part of the effective magnetic flux, B reg , only de¬ 
pends on the radius p and is spatially confined to the 


m = 0 — m = — 2 — m = —3 


1 

0.8 

§ 

P 0.6 

CjO 

to 

0.4 


0.2 
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FIG. 3: Magnon spectrum in the presence of a single skyrmion 
excitation as a function of k 2 /Q 2 measuring the strength of 
the magnetic field! 29 The magnon gap £ gap = £DMft 2 /Q 2 in¬ 
creases linearly with the field (black solid line). The field- 
polarized state becomes unstable at ft 2 r ~ 0.8Q 2 (dashed- 
dotted line) while the theory |l4| becomes locally unsta¬ 
ble at ftbimeron ~ 0.56Q 2 . Apart from the zero mode (not 
shown), there exist three subgap modes with angular momen¬ 
tum m — 0, —2, —3. 



skyrmion area. Its spatial distribution can be related 
with the help of Eq. ([9| to the topological charge den¬ 
sity p£ op of the skyrmion in addition to a term propor¬ 
tional to Q. While — p^ op is always positive, the lat¬ 
ter term can also be negative so that B reg as a func¬ 
tion of distance p even changes sign for lower values 
of ft 2 , see Fi g. |2| The spatial integral over the second 
term of Eq. (|24| ) however vanishes so that the total reg¬ 
ular flux f d 2 rB reg (p) = —Anh f d 2 rp^ op = Anh is quan¬ 
tized and det ermined by the topological charge of the 
skyrmion ! 25 * 49 * 


C. Magnon spectrum 

In order to solve Eq. ( [l4| ) for the magnon eigenval¬ 
ues and eigenfunctions, one uses the angular momentum 
basis T(r,£) = e~ l£t ^ h+irnx f] rn (p) with positive energy 
£ > 0. The angular momentum hm turns out to be a 
good quantum number and the wave equation ( p~4| ) re¬ 
duces to a radial eigenvalue problem for ffm{p) that can 
be solved with the help of the shooting method. 2 ^ In 
order to obtain positive expectation values of the Hamil¬ 
tonian, one has to look for eigenfunctions with a positive 
norm, 


dpprf m {p)T z rf m {p) > 0. (25) 

The resulting spectrum is shown in Fig. [3] as a func¬ 
tion of the parameter ft 2 /Q 2 that measures the strength 
of the magnetic field. The magnon continuum with the 
scattering states are confined to energies larger than the 






























5 


magnon gap £ gap oc ft 2 , which increases linearly with the 
field (black solid line). In the field range shown, there are 
three subgap states that correspond to bound magnon- 
skyrmion modes. While the breathing mode with angu¬ 
lar momentum m = 0 exists over the full field range, a 
quadrupolar mode with m = — 2 emerges for lower fields 
just before the field-polarized state becomes globally un¬ 
stable at ft 2 r ~ 0.8Q 2 (dashed-dotted line). The eigenen- 
ergy of the latter finally vanishes at ftb imeron ~ 0.56Q 2 , 
indicating a local instability of the theory with respect to 
quadrupolar deformations of the skyrmion, i.e. the for¬ 
mation of a bimeron . 50 Furthermore, a sextupolar mode 
with m = —3 only exists within the metastable regime. 
The corresponding eigenfunctions of these modes do not 
possess any nodes. 

Apart from the modes shown in Fig. [3j the spectrum 
of T-L also contains a zero mode with angular momentum 
m = — 1 given by 

-| / sin 0 _ Q r \ 

= vi ( sJ ¥+ e ')' (26) 

This zero mode is related to the translational invariance 
of the theory © that is explicitly broken by the skyrmion 
solution. The real and imaginary part of the amplitude of 
the eigenfunction ( [26] ) correspond to translations of the 
skyrmion within the two-dimensional plane. We have 
not yet found bound modes with a single or more nodes, 
which might however emerge for m = —1 at larger fields. 


III. HIGH-ENERGY SCATTERING OF 
MAGNONS 


The properties of the magnon scattering states for arbi¬ 
trary energies, e > £ gap , have been discussed in Ref. [29j 
In the present work, we elaborate on the scattering of 
magnons in the high-energy limit, 6 £ gap , which cor¬ 
responds to magnon wavevectors much larger than the 
inverse skyrmion radius, kr s 1. In this limit, the treat¬ 
ment of the scattering simplifies considerably allowing for 
a transparent discussion of characteristic features. 

In the high-energy limit the magnon-skyrmion inter¬ 
action is governed by the scattering vector potential 
a( r) = a x (p)x of Eq. (19) so that the scattering has a 
purely magnetic character. In particular, in this limit 
one can neglect the anomalous potential V x , and the BdG 
equation (14) reduces to a Schrodinger equation for the 


magnon wavefunction 


, fcn , fh 2 (—iV — a ) 2 


J m a g 


gap^ V ; - 


(27) 


Setting ^(r, £) = e 


— p -i£kt/h p imx 


e irnx r]m(p) with the dispersion 


> + 


h z k z 


£k ~ 5 §ap t- 2M ln<1& 

radial wave equation for 77 m (p) 


and wavevector k > 0 , one obtains the 




9 P \ , ( m-pa x (p )) 2 


- k 2 


rj m = 0 . (28) 



bQ 


FIG. 4: Scattering phase shift for high-energy magnons (30) 
as a function of impact parameter b for different values of 
ft 2 /Q 2 . The scattering is non-perturbative as the phase shift 
assumes values within the entire interval (—7r, 7r). 


For large distances pa x (p) —)> 1, which identifies the angu¬ 
lar momentum of the incoming wave to be L z = h(m— 1). 


A. Eikonal approximation 


As we are interested in the high-energy limit, we can 
treat this wave equation in the eikonal approximation. 
However, in order to make contact with Ref. f29l we first 
give the resulting phase shift within the WKB approxi¬ 
mation that is obtained by following Langei^^ 


jwkb = j°° _ k y p 

+ 7 ^\m-l\-kp 0 (29) 

where po is the classical turning point. The eikonal ap¬ 
proximation for the phase shift is then obtained by tak¬ 
ing the limit k 00 while keeping the impact parameter 
b = L z /(hk ) fixed, ^ KB yielding 



lab 


0) 


vV 2 - & 2 


dp = b 


f; 


lab 


(s\ b \) 




ds 


(30) 


where we used paf ah (p) = pa x (p) — 1, see Eq. ( [22] ) , and 
in the last equation we substituted s = p/\b\. This phase 
shift is odd with respect to 6, i.e. 800 (b) = — 8 oo(—b). Note 
that the scattering is non-perturbative even in the high- 
energy limit in the sense that the phase shift 800 (b) covers 
the entire interval (— 7 r, 7 r) as a function of 6, see Fig. B In 
particular, in the limit of small impact parameter 6 —0: 


800 (b)-tb[ ‘ 2 ^==- ds = —Trsgn(fr). (31) 

J i Vs 2 - 1 

For impact parameters larger than the skyrmion radius, 
b r s , the phase shift vanishes exponentially. 
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B. Differential cross section 


In the following, we consider a magnon scattering setup 
where an on-shell magnon plane wave with wavevector 
k = kx along the x-direction and amplitude A defined 
within the laboratory orthogonal frame, see Eq. (12), is 
impinging on the skyrmion, see also Fig. [ljb). At large 
distances this wavefunction assumes the asymptotic be¬ 
havior 


/ f e ikp \ 

V>iab(M) = Ae-^ h (^e ikr + fix)— j 


(35) 


where the scattering amplitude is given by 


FIG. 5: Classical deflection angle for scattering of high-energy 
magnons ( |33| as a function of impact parameter b for different 
values of kJQ 2 . In the high-energy limit, the scattering is in 
the forward direction with a deflection angle decreasing with 
increasing wavevector k as ©!J g (fr) ~ 1 /k. The inset focuses 
on the change of curvature at b = 0 for n 2 « 1.6Q 2 with the 
same units on the vertical axis. 


fix) = 


0 —iir/4 


V2nk m= 


E 


s i(m-l)x^ e «2(5 Tr 


-i). 


(36) 


m= — oo 


Note that the additional phase factor e~ %x arises from the 
gauge transformation (12). The differential cross section 
is then obtained by = |/(x)| 2 . 


The deflection angle in the eikonal approximation is 
given by the derivative of 5^ (6), 


1. High-energy limit of the scattering amplitude 


0 oo ib) = = l^ib) = 0£ g (6) - ^5(b). 

(32) 

The step of 600 (b) for head-on collisions, see Eq. 
leads to the delta function 5(b). The classical deflection 
function is given by the regular part, which reads 


= 


A r°° s_m^(s\b\i ds 

hk J 1 \J s 2 - 1 
- J_j ree (V^)d X , 


(33) 

(34) 


where in the last equation we substituted x = |6|\A 2 — 1 
and used that the integrand is an even function of x. It is 
determined by the regular part of the flux density, S reg , 
given in Eq. (24), integrated along a straight trajectory 
shifted from the x-axis by the impact parameter b. Its 
behavior as a function of b is shown in Fig. [5] for various 
values of n 2 /Q 2 . The deflection angle is always positive 
implying that, classically, the Lorentz force attributed to 
S reg always skew scatters the magnons to the right-hand 
side from the perspective of the incoming wave even for 
negative impact parameters, see Fig. 0 b )- Note that 
the deflection angle possesses a local minimum at b = 0 
for n 2 < 1.6Q 2 , that however gets filled and transitions 
into a maximum for larger values of n. This change of 
curvature at b = 0 is related to the change of curvature 
of the flux density B' r ' eg (p ) at the origin p = 0, see Fig. [2j 
that happens for a similar value of n. As the total flux 
of S reg is quantized, the deflection angle integrated over 
the impact parameter is just given by the universal value 
H 00 dbQ^{b) = ^/k. 


In the high-energy limit, we can replace the sum over 
angular momentum numbers by an integral over the im¬ 
pact parameter, b = (m — l)/fc, so that the scattering 
amplitude reads approximately 


p — i7r/4 r oo 

/ 00 ( X ) = -==fc/ dbe ibk x(e i 2 S °°W-1), (37) 

V27r/c J- oo 


with £oo(fr) defined in Eq. (30). The differential cross 
section in this limit, 


ddry 


dx 


= l/oo(x)| 2 = -^2 S ( k x/Q), (38) 


is then determined by the dimensionless function S', 
which is shown in Fig. [6] 

The support of the differential cross section is approx¬ 
imately limited by the extremal values of the classical 
deflection angle of Eq. (33) and Fig. [5] Note that the 
angle x is defined in a mathematically positive sense so 
that a positive 0 translates to a negative value of y. It 
is strongly asymmetric with respect to forward scatter¬ 
ing reflecting the skew scattering arising from the Lorentz 
force of the emerging magnetic field £> reg . 


2. Rainbow scattering and Airy approximation 

Moreover, the differential cross section exhibits oscil¬ 
lations. These can be attributed to an effect known as 
rainbow scattering. As the function 0^ s (6) is even in 6, 
there exist for a given classically allowed deflection angle 
0 always at least one pair ±b c \ of impact parameters that 
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FIG . 6: Differential cross section of high-energy magnons 
(38) for various values of k 2 /Q 2 . It is asymmetric with re¬ 
spect to x = 0 due to skew scattering, and the oscillations 
are attributed to rainbow scattering. The inset compares the 
curve for n 2 /Q 2 = 2 with the Airy approximation (|39| (green 


solid line) with the same units on the vertical axis; the ar¬ 
row indicates the position of the corresponding rainbow angle 
-k®z s (o)/Q. 


FIG. 7: Total scattering cross section of the skyrmion in the 
high-energy limit, Eq. (40), as a function of n 2 /Q 2 . It de¬ 
creases for increasing external magnetic field strength k 2 . 


effect.^ Finally, for smaller values of k 2 there also ex¬ 
ist two pairs of classical trajectories that interfere in the 
differential cross section. 


solve 0^ g (±6 c i) = 0. For a given angle 0 the magnons 
might, therefore, either pass the skyrmion on its right- or 
left-hand side; these classical trajectories interfere lead¬ 
ing to the oscillations in dcr/dx- 

First, consider values n 2 > 1.6Q 2 for which 0J^ S (6) 
possesses only a single maximum at b = 0. The maximum 
value 0^ S (O) is known as rainbow angle and for values 
of x close to — 0^ S (O), the interference effect of classical 
trajectories can be illustrated with the help of the Airy 
approximation for the scattering amplitude. For such 


values of x> the —1 in the integrand of Eq. (37) can be 


neglected as it only contributes to forward scattering. 
Expanding the exponent of the remaining integrand up 
to third order in b one then obtains 


/oo(x) 


Airy 


(39) 


— Z7T1 4 f* 7 

( rfc J db exp i6fc(x + 0^ s (O))+i-0"^ s (O)6 3 

— OO 


a/ 27t/c 

y/2nk e~ i7T / 4 
[fc|0 ,, ^(O)|/2] 1 /3 


Ai - 


fc(x + 9£ s (0)) A 

[k\e"z g mm l/3 ^ 


where in the last equation we identified the integral repre¬ 
sentation of the Airy function Ai using that 0"^ S (O) < 0. 

In the inset of Fig. [6j we compare the differential cross 
section at k? = 2 Q 2 with the Airy approximation result¬ 
ing from Eq. (39). The latter reproduces the exponential 


decrease for large angles x < — 0^ S (O) corresponding to 
the dark side and also the oscillations on the bright side, 
X > —0^ S (O), of the rainbow angle. It of course fails 
close to forward scattering and for positive angles x > 0 
where the classical deflection angle has lost its support. 

Close to n 2 ~ 1.6Q 2 even the derivative 0"^ S (O) van¬ 
ishes, see inset of Fig. [5j giving rise to a cubic rainbow 


C. Total and transport scattering cross section 

We continue with a discussion of the total, cr to t = 
fl^dxdcr/dx, and the transport scattering cross section 
defined in Eq. In order to determine their high- 

energy limit, one first expresses dcr/dx = |/(x)| 2 in terms 
of the exact representation ( [36] ) for the scattering ampli¬ 
tude /(x) and evaluates the integral over x- Afterwards 
one takes the high-energy limit k -A oo with keeping the 
impact parameter b = (m — l)/k fixed. 

The total scattering cross section of the skyrmion then 
reduces to 



(40) 


It saturates to a finite value in the high-energy limit, 
and its dependence on k is shown in Fig. [7] It decreases 
with increasing n and thus decreasing skyrmion radius 
r s as expected. One might expect that a^ t ~ r s which 
however only holds approximately. 

Using that S^^b) is an odd function of 6, we obtain 
for the transport scattering cross section a l(^) in the 
high-energy limit 



<5^ (6) (sin (5oo (6)) 2 = 

_ sin(25oo) 1 0 _ 4it 

4 - —7T k 


(41) 

(42) 


In the last line, we further used the boundary values of 
the function 6 ^(b). It vanishes cr^(e) ~ 1/fc, but with a 
universal prefactor that is independent of k. 

Finally, for the ongitudinal transport scattering cross 





























FIG. 8: The longitudinal transport scattering cross section, 
Eq. (44), vanishes as a y° ~ 1/k 2 in the high-energy limit. 
The panel shows the /^-dependence of the prefactor. 


section we obtain for kr s 1 


a 


oo 

II 


(4 


^2 j db ( 2 (0 2 ( sin<5 oo) 2 
0 


- 5'^ sin <*«, 



(43) 


After integrating by parts this simplifies to 


are given by 


F j_ = ka^(e)(z x J £ ) = 47r(i x J £ ) = —G x iT e fr , (46) 
F„ = ka^(e)J £ = F db(eZ g (b)) 2 v eS , (47) 

J — C 


87T 


where we used Eqs. (41) and (44) as well as G — —\G\z. 
They are indeed of the form given in Eq. ©• The 
transversal momentum-transfer force, Fj_, is universal, 
and F\\ is determined by the (3 parameter of Eq. © af¬ 
ter identifying 0(6) with the classical deflection angle 

©sw- 

Is there an intuitive classical interpretation of these 
momentum-transfer forces? From the classical limit of 
the Schrodinger equation ([27]) follows the equation of 
motion for the coordinate r[t) of a classical magnon 
particl^l 


M mag r = fx (zB reg (\r\)) 


(48) 


with the regular part of the effective magnetic flux dis¬ 


tribution B reg of Eq. (24). Note that we have chosen in 


Eq. (27) the charge to be +1. Consider the change of 


momentum, of this magnon particle after scattering 
off the static skyrmion by integrating the left-hand side 
of Eq. (48), 


4 r°° r°° i 

= ¥ J db(6’ oc (b)f= y db-(0S?(6)) a 


(44) 


It is given by the square of the classical deflection angle 
(33) integrated over the impact parameter b. It vanishes 
as a j^° ^ 1/k 2 in the high-energy limit with a prefactor 
whose k dependence is shown in Fig. [8] On dimensional 
grounds one might expect k 2 a^ ~ l/r s , which again 
only holds approximately. 


D. Magnon pressure in the high-energy limit 


We have shown in Ref. [29] by considering the energy- 
momentum tensor of the field theory that the monochro¬ 
matic plane wave of (35) with wavevector k = kx leads 


to a momentum-transfer force in the Thiele equation of 
motion of the form given in Eq. ©> with the magnon 
current 


Je 


= x\A \ 2 


moh Hk 
g/i B M mag 



(45) 


In the second equation, we have introduced the effective 
velocity v e q = x\A\ 2 and \G\ = AttuioH/( gpB) with 
the purpose of comparing with Eq. ©>• 

This momentum transfer is illustrated in Fig. [9] In the 
high-energy limit, the transversal and longitudinal forces 


5p(b) = / dtM mag r(t ) M mag (f(oo) - f(-oo)) 

J — oo 

= „( cos0(6)-l\ 

^ f —sin 0(b) J ' ' 

In the last equation, we have exploited that at large dis¬ 
tances the magnitude of momentum M mag |r(±oc)| = p 
remains unchanged due to energy conservation, while the 
orientation of velocity is determined by the scattering 
angle 0(6), see Fig. [ljb), that depends on the impact 
parameter 6 of the trajectory. 

This momentum 5p(b) is transferred to the skyrmion. 
The momentum-transfer force on the skyrmion due to a 
current of classical magnon particles along x with density 
m oK qI^b) and velocity v e s = \v e ^\ is then given by 

F = ( p 11 ) = ~ v eS-^~ [ db5p(b), (50) 
J gg B i-oo 


with F\ |/_l = \Fu /_l|. In the high-energy limit, the scatter¬ 
ing is in forward direction so that we can expand Eq. (49) 
in the deflection angle 0(6) and the force becomes with 
p = hk 


r = v e ff - hk 

9hB 


F 

J —c 


db 


mb)) 2 

0 ( 6 ) 


(51) 


Finally using that the integral dbO(b) = Air/k is 
quantized in the high-energy limit, that we already know 
from the discussion in the context of Eq. ( [33] ) , we recover 
Eqs. (46) and (47). 



















9 



FIG. 9: An incoming monochromatic magnon current J £ leads 
to a momentum-transfer force F that is determined by the 
transport scattering cross sections, see Eq. The image 
shows the magnon wavefunction in the WKB approximation 
with the skyrmion being represented by the circle with radius 
r s [29]. For high-energy magnons with wavevector kr s 1, 
the transversal force dominates, F\\/F± ~ 1/A;, resulting in a 
skyrmion motion dtR approximately antiparallel to J £ with a 
small skyrmion Hall angle <f> ~ 1/k. 


For the understanding of the universality of F l* it is 
also instructive to consider alternatively the right-hand 
side of the classical equations of motion (48). By inte¬ 


grating the right-hand side, one obtains for the transver¬ 
sal momentum change 

/ OO f'OO 

dt(-x)B ies (\r\) « - / <M? reg (\/& 2 + a; 2 ). 

-oo J —OO 

(52) 


In the last equation we employed the high-energy ap¬ 
proximation by straightening the magnon trajectory. It 
follows then for the transversal force 

/ OO fOO 

db / dxt3 reg (\/b 2 + x 2 ) (53) 

-oo J — oo 

= v eS ^-4nh, (54) 

9 \*b 

where its universality is now directly related to the quan¬ 
tized total flux Of S reg . 


IV. SUMMARY 

The scattering of high-energy magnons with wave vec¬ 
tors kr s ^ 1 off a magnetic skyrmion of linear size r s 
is governed by a vector scattering potential. The associ¬ 
ated effective magnetic field is related to the topological 
charge density of the skyrmion and is exponentially con¬ 
fined to the skyrmion area. The total flux is determined 
by the topological skyrmion number and is quantized. 


When a magnon traverses the skyrmion, classically 
speaking, it experiences the resulting Lorentz force and is 
deflected to a preferred direction determined by the sign 
of the emergent magnetic flux. This results in skew scat¬ 
tering with a differential cross section that is asymmetric 
with respect to forward scattering, see Fig. [6] As the 
flux distribution is rotationally symmetric, the classical 
deflection angle 0(6) as a function of the impact param¬ 
eter 6 is even in the high-energy limit, 0(6) = @(—6). As 
a consequence, for a given deflection angle 0 there exist 
corresponding classical trajectories with positive as well 
as negative 6, i.e., that pass the skyrmion on the left-hand 
as well as on the right-hand side. These trajectories in¬ 
terfere which leads to oscillations in the differential cross 
section, an effect known as rainbow scattering. 

Magnons hitting the skyrmion also transfer momentum 
giving rise to a force in the Thiele equation of motion, see 
Eq. (©• In the high-energy limit, this force can be inter¬ 
preted classically and assumes the form of Eq. <§• While 
the transversal momentum-transfer force, F± is universal 
and determined by the total emergent magnetic flux, the 
longitudinal momentum-transfer force, Fy is obtained by 
integrating (0(6)) 2 over the impact parameter 6 leading 
to the parameter /3 £ of Eq. <©• Since for large energies 
the classical deflection angle is small, 0(6) ^ 1/fc, the 
momentum transfer is mainly transversal, F\\/F± ~ 1/k. 

This leads to a skyrmion motion dtR approximately an¬ 
tiparallel to the magnon current J £ with a small skyrmion 
Hall angle 4> = /3 e /\G\ defined in Fig. [ 9 J 


4 > = 


1 szjmfdb 

2 fZmdb 



(55) 


where the integral J'^° Q(b)db = 4ir/k is universal in the 
high-energy limit. Interestingly, the Hall angle 4> at high 
energies increases with increasing ft, which is shown in 
Fig. |8j identifying 4> = ka^ (s) / Fk . 

While the skyrmion Hall angle 4> is small at high en¬ 
ergies kr s 1, we note that it increases with decreasing 
energy and assumes the maximum value 29 <f> = 7r/2 in 
the low-energy limit kr s <C 1 where 5-wave scattering 
prevails and Eq. § ceases to be valid. 
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